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^ \ Abstract 

Q_i' We report preliminary results on the chiral and deconfinement aspects of the 

^ . QCD transition at finite temperature using the Highly Improved Staggered Quark 

(HISQ) action on lattices with temporal extent of N T = 6 and 8. The chiral aspects 
■ of the transition are studied in terms of quark condensates and the disconnected 

chiral susceptibility. We study the deconfinement transition in terms of the strange 
quark number susceptibility and the renormalized Polyakov loop. We made con- 
tinuum estimates for some quantities and find reasonably good agreement between 
O ■ our results and the recent continuum extrapolated results obtained with the stout 

D \ staggered quark action. 



Introduction 



Improved staggered fermion formulations are widely used to study QCD at non-zero 
temperatures and densities, see e.g. Ref. jU |2] for recent reviews, for, at least, two 
reasons: they preserve a part of the chiral symmetry of continuum QCD which allows 
one to study the chiral aspects of the finite temperature transition, and are relatively 
inexpensive to simulate numerically because due to the absence of an additive mass 
renormalization the Dirac operator is bounded from below. However, lattice artifacts 
related to taste symmetry breaking turned out to be numerically large. To reduce the 
taste violations smeared links, i.e. weighted averages of different paths on the lattice 
that connect neighboring points, are used in the staggered Dirac operator and several 
improved staggered formulations, like p4, asqtad, stout and HISQ differ in the choice of 
the smeared gauge links. The ones in the p4 and asqtad actions are linear combinations 
of single links and different staples [H H] and therefore are not elements of the SU(3) 
group. 

It is known that projecting the smeared gauge fields onto the SU(3) group greatly 
improves the taste symmetry [5]. The stout action [6] and the HISQ action implement 
the projection of the smeared gauge field onto the SU(3) (or simply U(3)) group and 
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thus achieve better taste symmetry at a given lattice spacing. For studying QCD at high 
temperature it is important to use discretization schemes which improve the quark dis- 
persion relation, thus eliminating the tree level 0(a 2 ) lattice artifacts in thermodynamic 
quantities. The p4 and asqtad actions implement this improvement by introducing 3-link 
terms in the staggered Dirac operator. 

In this paper we report preliminary results on the chiral and deconfinement transition 
in QCD at non-zero temperature obtained with the HISQ action which combines the 
removal of tree level 0(a 2 ) lattice artifacts with the addition of projected smeared links 
that greatly improve the taste symmetry. We also compare our results to the continuum 
extrapolated results obtained with the stout action [7]. 



1 Action and run parameters 

The Highly Improved Staggered Quark (HISQ) action developed by the HPQCD/UKQCD 
collaboration [8] reduces taste symmetry breaking and decreases the splitting between 
different pion tastes by a factor of about three compared to the asqtad action. The net 
result, as recent scaling studies show [HI ID], is that a HISQ ensemble at lattice spacing a 
has scaling violations comparable to ones in an asqtad ensemble at lattice spacing 2a/3. 

In this study we used the HISQ action in the fermion sector and the tree-level 
Symanzik improved gauge action. The strange quark mass m s was set to its physi- 



cal value adjusting the quantity \J2rn\ — m 2 = m Vsg ~ y/2Bm s to the physical value 
686.57 MeV. We used two values of the light quark mass: mi = 0.05m s and mi = 0.20m s . 
These correspond to the lightest pion mass of about 160 MeV and 320 MeV respectively 
[TTj . Calculations have been performed on 24 3 x 6 and 32 3 x 8 lattices for the smaller 
quark mass, while for the larger quark mass we used 16 3 x 6 lattices. At several values of 
the lattice spacings zero temperature simulations have been performed on 32 4 lattices for 
mi = 0.05m s and on 16 3 x 32 lattices for mi = 0.2m s . The molecular dynamics (MD) tra- 
jectories have length of 1 time unit (TU) and the measurements were performed every 5 
TUs at zero and 10 TUs at finite temperature. However, for the few smallest beta values 
the trajectories had length of 1/2 or 1/3 TU. In this case measurements have been per- 
formed after each 20 or 30 trajectories (which corresponds, again, to 10 TUs). Typically, 
at least 300 TUs were discarded for equilibration at the beginning of the simulations. 

The lattice spacing has been determined by measuring the static quark anti-quark 
potential. As in previous studies by the MILC collaboration the static potential was 
calculated fixing to Coulomb gauge and considering temporal Wilson lines of different 
extent. Forming the ratio of these correlators and fitting them to a constant plus linear 
Ansatz we extracted the static potential V(r). We have calculated the Sommer scale r 
defined as 



2 dV 
dr 



1.65. (1) 
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The potential has been normalized to the string potential 



7T 

V string = ~~^r~ ° T ' ^ 



2 




110 120 130 140 150 160 170 180 190 200 120 140 160 180 200 220 



Figure 1: The subtracted chiral condensate and A/* compared to the continuum estimate 
for the stout action |7J. The solid black line is our continuum estimate for the HISQ 
action. 

at r = 1.5r or equivalently to the value 0.91/r at r = r . The additive constant 
determined by this normalization is used to calculate the renormalization constant for 
the Polyakov loop as will be discussed later. To convert from lattice units to physical 
units we use the value = 0.469 fm as determined in [12]. 

2 Chiral and deconfinement transition in QCD 
2.1 The chiral transition 

In the limit of zero light quark masses QCD has a chiral symmetry and the finite tem- 
perature transition is a true phase transition. The order parameter for this transition is 
the light chiral condensate (^pip)i- However, even at finite values of the quark mass the 
chiral condensate will show a rapid change in the transition region indicating an effective 
restoration of chiral symmetry. The fluctuations of the chiral condensate, also called the 
disconnected chiral susceptibility, will have a peak at the transition temperature. Recent 
studies with the p4 action suggest that for the physical strange quark mass and two 
light quarks the chiral transition is of the second order for vanishing light quark mass, 
belonging to the O(N) universality class [13]. Thus, the universal properties of the chiral 
transition in the limit of zero light quark masses govern the transition for sufficiently 
small but non-zero light quark masses [13]. The corrections to scaling turned out to 
be small for mi = m s /20 and thus the temperature variation of the chiral condensate 
is to a large extent determined by the singular part of the free energy density which is 
universal [13]. This allows to determine the chiral transition temperature for non-zero 
quark masses. The temperature derivative of the chiral condensate and disconnected 
chiral susceptibility will diverge in the limit m; — > 0. Since the chiral condensate has an 
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Figure 2: The disconnected chiral susceptibility calculated with the HISQ action on 
N T = 6 and 8 lattices for mi = 0.05m s and mi = 0.20m s . 

additive ultraviolet renormalization we consider the subtracted chiral condensate [23] 

Here the subscript / and s refer to light and strange chiral condensates, while subscript 
and t to the expectation value at zero and non-vanishing temperature. Another pos- 
sibility to get rid of the additive renormalization in the chiral condensate is to consider 
the quantity 

Af (T) = -m„fft$V),, T - (^>*,o) (4) 

which up to a constant multiplicative factor is identical to the quantity introduced in Ref. 
[7] and was called the renormalized chiral condensate. In Fig. [Tjwe show our numerical 
results for A^ s and Af and compare them to the continuum extrapolated stout results 
[7]. To facilitate this comparison for mi = 0.05m s we performed a combined polynomial 
fit of the N T = 6 and N T = 8 results allowing for lattice spacing dependent coefficients 
in this fit. This allows us to give an estimate of these quantities in the continuum limit 
shown as the solid black lines. The continuum estimates for HISQ are different from 
the continuum extrapolated stout results. This is presumably due to the fact that the 
quark masses used in the stout calculations are smaller, namely mi = m s /27.3 [H]. The 
transition temperature is decreasing with decreasing the quark mass. Also the chiral 
condensate at fixed temperature decreases with decreasing the quark mass. Thus at a 
qualitative level this discrepancy could be understood. Future calculations on iV T = 12 
lattices and smaller quark masses will be needed to clarify this issue quantitatively. 

We also calculated the disconnected chiral susceptibility. In Fig. [2] we show the light 
quark disconnected chiral susceptibility. The peak position in this quantity also defines 
the chiral transition temperature. To estimate the peak position in this quantity we have 
performed fits of the lattice data using different functional forms. We also varied the fit 
intervals. Our analysis gives 

N T = 6 : T c = 168(4)(5)MeV (0.05m s ); T c = 185(4)(5)MeV (0.20m s ) (5) 
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N T = 8 : T c = 165(4)(5)MeV (0.05m s ), 



(6) 



where the first error reflects the uncertainty in the determination of the peak position 
estimated using different fit forms and varying the fit range. The second error is due to 
the scale determination. The peak positions in the disconnected chiral susceptibilities are 
significantly lower than those obtained with p4 and asqtad actions for N T < 8 [HI [17] , 
though they are compatible with the recent estimates for asqtad on iV T = 12 lattices 
|llj . Unfortunately no published stout data are available for the disconnected chiral 
susceptibility and direct comparison between HISQ and stout results is not possible here. 
Therefore, in Ref. [TTJ we compared HISQ and stout data for the renormalized chiral 
susceptibility introduced in [15]. Some small discrepancies between the stout and HISQ 
results have been found there which again are probably due to slightly different quark 
masses. Note that within errors of the calculations the peak positions in the renormalized 
chiral susceptibility are consistent between HISQ and stout [TTj . 

2.2 Deconfinement transition 

By the deconfinement transition we mean liberation of many degrees of freedom, which 
also could be understood as a transition from hadronic degrees of freedom to partonic 
ones, and the onset of color screening. The aspects of the deconfinement transition related 
to color screening are studied in terms of the Polyakov loop and Polyakov loop correlators. 
The Polyakov loop needs to be renormalized and after proper renormalization it is related 
to the free energy of a static quark anti-quark pair at infinite separation F^iT) [18], 



The renormalized Polyakov loop can be obtained from the bare Polyakov loop as 



Here the multiplicative renormalization constant z((3) is related to the additive normal- 
ization of the potential c(/3) as z((3) = exp(— c(/3)/2) discussed in the previous section. 
The renormalized Polyakov loop has been calculated for pure gauge theory [TSJ HH] , 3- 
flavor QCD [20J as well as for 2-flavor QCD [2TJ. More recently it has been calculated 
for 2+1 flavor QCD with physical strange quark mass and light quark masses close to 
the physical values [7J [231 [22 [251 HS1 EH] ■ In Fig. [3] we show our results for the renor- 
malized Polyakov loop calculated with the HISQ action using N T = 6 and 8 lattices and 
mi = 0.05m s and compare them to the continuum extrapolated results obtained with 
the stout action |7J. The decrease of F^T), and thus the increase in the Polyakov loop 
could be related to the onset of screening at high temperatures (e.g. see discussion in 
Ref. [22]). On the other hand, in the low-temperature region the increase of L ren is 
related to the fact that there are many static-light meson states that can contribute to 
the static quark free energy close to the transition temperature, while far away from the 
transition temperature it is determined by the binding energy of the lowest static-light 



L ren (T) = exp(-F 00 (T)/(2T)). 



(7) 
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Figure 3: The renormalized Polyakov loop (left) and the strange quark number suscep- 
tibility (right) calculated with the HISQ action for mi = 0.05m s and compared to the 
continuum extrapolated stout results. The lines correspond to continuum extrapolation. 



mesons. The strange quark number susceptibility is defined as the second derivative of 
the pressure with respect to the quark chemical potential 

_ &p(T, fl8 ) i 

Xs g ^ U=o- W 

It describes strangeness fluctuations at zero strange quark chemical potential. At 
low temperatures strangeness is carried by strange hadrons, which are heavy compared 
to the temperature. As a result strangeness fluctuations are suppressed in the low- 
temperature region. At high temperatures, on the other hand, strangeness is carried by 
light quarks and strangeness fluctuations are close to the value given by an ideal quark 
gas. Deconfinement will manifest itself as a rapid increase in the strangeness fluctuations 
in some temperature interval, reflecting the change in the relevant degrees of freedom 
from hadronic to partonic. Therefore, strangeness fluctuations are used as a probe of 
the deconfinement transition [231 Ell [13 EH] • In Fig. [3] we show the strange 
quark number susceptibility calculated for the HISQ action on N T = 6 and 8 lattices for 
mi = 0.05m s . 

We performed a combined polynomial fit of our lattice results where the coefficient of 
the polynomial had a 2 correction. This allows to estimate strangeness fluctuations in the 
continuum limit. We compared our results to those obtained with the stout action and 
extrapolated to the continuum. As one can see there is reasonable agreement between 
HISQ and stout results in the continuum limit. The main difference is the inflection 
point in the stout calculations at temperatures of about 165 MeV which is not visible 
in the HISQ results. Note, that strangeness fluctuations approach the continuum limit 
from below. This is expected due to the lattice spacing dependence of the hadron masses 
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3 Conclusions 



We studied the deconfinement and chiral transition in QCD at non-zero temperature 
using the HISQ action and compared our results to recent continuum extrapolated results 
obtained with the stout action. We find reasonable agreement in the continuum limit 
between the results obtained with different actions for the renormalized Polyakov loop 
and the strange quark number susceptibility. We also calculated the chiral condensate 
and disconnected chiral susceptibility for the HISQ action. For these quantities a direct 
comparison with stout results is more difficult due the slightly different light quark masses 
used in the two calculations. Nonetheless, current HISQ results indicate a chiral transition 
temperature that is significantly lower than the previous estimates obtained with the 
asqtad and p4 actions using lattices with the temporal extent N T < 8. It would also be 
interesting to compute the equation of state with the HISQ action and compare it to the 
very recent results obtained with the stout action [27]. 



Acknowledgements 

This work has been supported in part by contracts DE-AC02-98CH10886 and DE-FC02- 
06ER-41439 with the U.S. Department of Energy and contract 0555397 with the National 
Science Foundation. The numerical calculations have been performed using the USQCD 
resources at Fermilab as well as the BlueGene/L at the New York Center for Computa- 
tional Sciences (NYCCS). 

References 

[1] C. DeTar and U. M. Heller, Eur. Phys. J. A 41, 405 (2009) 

[2] P. Petreczky, Nucl. Phys. A 830, 11C (2009); P. Petreczky, Nucl. Phys. Proc. Suppl. 
140, 78 (2005) 

[3] F. Karsch, E. Laermann and A. Peikert, Nucl. Phys. B 605, 579 (2001) 

[4] K. Orginos, D. Toussaint and R. L. Sugar [MILC Collaboration], Phys. Rev. D 60, 
054503 (1999) 

[5] A. Hasenfratz, Nucl. Phys. Proc. Suppl. 119, 131 (2003) 
[6] Y. Aoki, et. al, JHEP 0601, 089 (2006) 

[7] S. Borsanyi et. al, [Wuppertal-Budapest Collaboration]. [arXiv:1005.350 8 [hep-lat]. 

[8] E. Follana et al. [HPQCD collaboration and UKQCD collaboration], Phys. Rev. D 
75 (2007) 054502 

[9] A. Bazavov et al. [MILC Collaboration], PoS LAT2009, 123 (2009) 
[10] A. Bazavov et al. [MILC Collaboration], lar2£iyil004 i Q342J [hep-lat]. 



7 



[11] A. Bazavov and P. Petreczky, J. Phys. Conf. Ser. 230, 012014 (2010) 

[12] A. Gray et al, Phys. Rev. D 72, 094507 (2005) 

[13] S. Ejiri et al, Phys. Rev. D 80, 094505 (2009) 

[14] Y. Aoki et al, JHEP 0906, 088 (2009) 

[15] Y. Aoki et al, Phys. Lett. B 643, 46 (2006) 

[16] M. Cheng et al, Phys. Rev. D 74, 054507 (2006) 

[17] Rajan Gupta. larXiv:0912l374l 

[18] O. Kaczmarek, et. al, Phys. Lett. B 543, 41 (2002) 

[19] S. Digal, S. Fortunato and P. Petreczky, Phys. Rev. D 68, 034008 (2003) 

[20] P. Petreczky and K. Petrov, Phys. Rev. D 70, 054503 (2004) 

[21] O. Kaczmarek and F. Zantow, Phys. Rev. D 71, 114510 (2005) 

[22] P. Petreczky, Eur. Phys. J. C 43, 51 (2005) 

[23] M. Cheng et al, Phys. Rev. D 77, 014511 (2008) 

[24] A. Bazavov et al, Phys. Rev. D 80, 014504 (2009) 

[25] M. Cheng et al, Phys. Rev. D 81, 054504 (2010) 

[26] P. Huovinen and P. Petreczky, Nucl. Phys. A 837, 26 (2010) 

[27] S. Borsanyi et al, [Wuppertal-Budapest Collaboration], larXiv:1007.2580l [hep-lat]. 



8 



